Vicinal Surfaces, Fractional Statistics and Universality 



On 
ON 
Os 

as: 

< 
o 



43 

o 



CO 

c3 



I 

C 

O 

o 



(N 
> 

m 
m 

o 

OS 
0\ 



13 

c 

o 
o 



X 
S3 



Somendra M. Bhattacharjeefcl 

Institute of Physics, Bhubaneswar 751 005, India 

Sutapa MukherjiEl 

Department of Physics and Meteorology, Indian Institute of Technology, Kharagpur 721 302, India 

(February 1, 2008) 

We propose that the phases of all vicinal surfaces can be characterized by four fixed lines, in 
the renormalization group sense, in a three-dimensional space of coupling constants. The observed 
configurations of several Si surfaces are consistent with this picture. One of these fixed lines also 
describes one-dimensional quantum particles with fractional exclusion statistics. The featureless 
steps of a vicinal surface can therefore be thought of as a realization of fractional-statistics particles, 
possibly with additional short-range interactions. 
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A crystal surface cut at a small angle to a symmetry 
direction is called a vicinal surface Q. Such a miscut 
surface consists of terraces of the symmetry plane, sep- 
arated by monatomic steps running across the sample 
in a preferred direction dictated by the cut. The den- 
sity of the steps (number of steps per unit transverse 
length), p, is related to the miscut angle (p ~ tan 6*). 
In thermal equilibrium, the configuration of the surface 
is determined by the steps and their interactions. Vicinal 
surfaces constitute a special class of objects in low dimen- 
sional statistical mechanics and are useful as substrates 
in many technological and experimental situations. An 
understanding of the phases and phase transitions of vic- 
inal surfaces is therefore important. 

Various types of behaviors of vicinal surfaces are 
known. For example, it is known that vicinal Si(lll) 
surfaces have single steps or triplets j|] but never pairs, 
while pairs of steps are seen on Si(001) jj). An un- 
usual phenomenon occurs in the case of Si(l 13) where a 
uniform-step-density phase, on cooling, phase-separates 
into a flat surface and a phase with a large angle (i.e., 
a high density of steps) |@|§. In other words, a small- 
angle vicinal Si(113) surface becomes thermodynamically 
unstable. The coexistence curve, on the temperature- 
versus-miscut-angle orientational phase diagram, ends at 
a tricritical point with the shape given by an expo- 
nent (3 = 0.42-0.54. A variant of the phase diagram has 
also been reported in Ref . [4b] . Is there a natural way of 
characterizing this zoo of vicinal surfaces? 

In this paper our aim is to develop a long-distance, uni- 
versal behavior for vicinal surfaces, where microscopic 
details like the specific material, the lattice structure, 
surface reconstruction if any, etc., do not play a direct 
role. This is best done in a continuum approach |^| pj. 
The steps are treated as fluctuating elastic strings (di- 
rected polymers), and all the effects of the surface go 
into the elastic energies of the steps and their effective 
interactions. Therefore, our approach is valid for any 



vicinal surface that can be characterized by featureless, 
wandering steps. 

A major interaction of the steps is known from the- 
ory ancl experiments pi to be | Vi{z) — rj(z) |~ 2 
where z is the special direction of the steps and rj is the 
transverse d-dimensional coordinate of a point at z along 
the length of the ith step. (Our interest is, of course, at 
d = 1.) This long-range interaction is generally repul- 
sive, originating from the elastic effects of the terraces, 
although, for metals, dipoles or quenched impurities on 
the steps can produce an attractive r~ 2 potential A 
theoretical explanation H of the observed tricriticality 
H requires an additional attractive short-range interac- 
tion. The continuum description of Ref. (7j is justified 
a posteriori by the existence of a renormalization group 
fixed point with a diverging length-scale (see also Refs. 
§0). 

In an effective Hamiltonian for the steps, the inter- 
action need not be restricted to just pairwise-additivc 
potentials, and, in a renormalization-group (RG) ap- 
proach, one should include allowed marginal operators. 
We keep the three-body short-range interaction since it is 
marginal at d = 1 (see below) |l2|,[7j . All m-step interac- 
tions are irrelevant for m > 3 at d = 1. The Hamiltonian 



v-25\{v i:j (z)) + 



H = H / *y*?w+E / dz 

+v 3 ^2 / dz8A,(Tij(z))6A(rjk(z)). 



(1) 



where Yi(z) = dvijdz, K{ is the elastic constant of the 
i th step, each of length N(— > oo), and g, V2, v 3 are re- 
spectively the long-range two-body, the short-range two- 
body and the short-range three-body coupling constants. 
The short-range interactions are taken as contact inter- 
actions, and power counting shows that this is sufficient 
p3[ . There is a short-distance cut-off (reminiscent of the 
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lattice) so that, in Fourier modes, <5a(q) = 1, for | q |< A, 
and 0, otherwise. For the tricritical point, one takes [QJ|] 
V2 = V20 (T — T t ) but, for generality, we consider both 
positive and negative values of all the three parameters 
(g, V2,Va). The partition function is then given by the 
summation of the Boltzmann factor over all possible con- 
figurations with free end-points. Two special cases of Eq. 
(0) ( v 3 — 0; an d g = 0) were considered in Ref. in 
connection with the tricritical point seen in Si(113). The 
physical picture used is that of phase separating polymers 
with the critical point of the coexistence curve coincid- 
ing with the binding-unbinding critical point of two steps 
(zero density). It was shown in Ref. that (3 = 0.5 
occurs for g — 3/4. 

The resemblance of Eq. (|l|) with the Calogero- 
Sutherland model @,^5| is apparent and we shall discuss 
some issues related to this. The Hamiltonian (U = 1) de- 
scribes a set of interacting quantum particles which can 
be chosen as bosons. The thermodynamic properties of 
the quantum system can be obtained from Eq. (|l]) with 
N as the inverse temperature. In the limit of infinite 
lengths of the steps, the properties of the vicinal surface 
are given by the ground state of the corresponding quan- 
tum problem ||. 

The main contents of the paper are the following: (a) 
We show that for, the v$ = case, at g = 3/4 the 
continuous transition goes over to a first-order transi- 
tion J]l5]jl^]. This requires a study of reunion exponents 
Jl7| , [l8| ] for the steps with the long-range interaction. The 
change in the order of the transition is then connected 
to the observed behavior of Si(113). (b) We obtain a 
fixed-point (or, rather, a fixed-line) description for the 
above Hamiltonian, and, from the nature of the fixed 
points, we argue that the observed features of various Si 
vicinal surfaces are described by these fixed points. We 
then conjecture that all vicinal surfaces are described by 
these fixed lines (and the flows) in a three-dimensional 
parameter space, (c) It is quite common to use a quan- 
tum description where the steps are treated as fermions, 
but we conclude, from the equivalence with the Calogero- 
Sutherland model on the ^3 = fixed line, that the steps 
should rather be treated as one- dimensional fractional- 
exclusion- statistics particles with only short-range inter- 
actions. 

So far as the RG flows of the various parameters are 
concerned, one may just restrict oneself to two or three 
steps ( "vertex functions" ) . For a many-step system, 
one also needs the chemical potential /Lt, but its flow is 
determined by dimensional analysis |I^,[l5|,|l9| . In the 
momentum- shell renormalization-group approach, short- 
wavelength fluctuations are integrated out in a thin shell 
(A/6, A) in Fourier space, and the effect is absorbed by 
redefining the parameters. The system is then rescaled 
to its original state, whereby the cutoff goes back to A. 
(We choose A = 1.) The change of the parameters with 
distance in the long-distance limit is then recovered from 



the RG flow-equations. The spatial rescaling factor is 
anisotropic, r — > 6r, z — > b^z, and the choice £ = 2 
keeps the elastic constant invariant For simplicity, 
we choose k, = 1, for all i. 

Since the RG transformation is analytic, the long- 
range term does not get renormalized. That there is 
no renormalization of g turns out to be rather natural 
when we discuss the quantum problem. We just quote 
the well-known recursion relations in terms of the dimen- 
sionless variable §112= L 2 - d K d v 2 +g/(d-2), where 
K d = 2ir d ' 2 /\r(d/2)(2ir) d ] and L is an arbitrary length- 
scale in the transverse direction (b = 1 + SL/L). Restrict- 
ing ourselves to d = 1, the flow equations are |l5|Jl9|,[l^| 



dui n dq da , , 

L— =u 2 -4 +9 , LJL=0, and L-Jjr — 2/x. (2) 

The flow equation for u 2 has two fixed points u* u — 
{1 ± [1 + 4g] 1 / 2 }/2, s,u denoting respectively the sta- 
ble and unstable fixed points (Fig. la). The unstable 
fixed point describes the binding-unbinding transition of 
two steps. The critical point has a length-scale exponent 
p± = l/\/l + 4<7 . For a many-step system this unstable 
fixed point corresponds to the tricritical point mentioned 
earlier, provided the fixed point describes a critical be- 
havior. Since the density p vanishes at this tricritical 
point, the assumption of a single length scale then tells 
us that p _1 should diverge as the length scale. Therefore, 
(3 = v±. For (3 = 1/2, as for Si(113), one would require 
g = 3/4 as obtained in Ref. 0. 

We now show that g = 3/4 is indeed a special point. 
Let us consider the free energy |l^] of two steps. Config- 
urations of the steps are of the type shown in Fig. 2. The 
bubble-like contributions (region B in Fig. 2a) represent 
the steps in the high temperature phase described by the 
stable fixed point. It is shown in Ref. jl7| (we, there- 
fore, skip the details) that the nature of the singularity 
of the free energy is determined by the decay of the parti- 
tion function of the bubble of length N, Z R {N) ~ N~^, 
where tp is the reunion exponent p8[ . If 1 < ip < 2, then 
the transition is continuous, while for 2 < ip it is first or- 
der, but with a weak singularity and a diverging length 
scale @. 

The reunion exponent ip determines the behavior of the 
partition function of two steps starting at the origin and 
reuniting anywhere in space fl"^ , |l8|| . We need the expo- 
nent in the high-temperature phase characterized by the 
stable fixed point. The zeroth- and the first-order dia- 
grams for the reunion partition function Z R are shown in 
Figs. ^|b,c. These contributions require renormalization 
of the partition function itself, over and above the renor- 
malization of the parameters. This extra renormalization 
yields a nontrivial ?/> different from the Gaussian value 
ip = d/2 in d dimensions. In the limit N — > 00, the loop 
contributions lead to the recursion relation for Zr(N) as 
L dZ R /dL = -2(d/2 + u 2 )Z R , so that at d = 1 at the 
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stable fixed point u\ = u s , we get ip = 1 + (1 + 4g) 1 / 2 /2. 
This is an exact result from which we recover the vicious 
walker exponent JlTj if) = 3/2 for g = 0. (The details 
and generalizations to arbitrary number of steps will be 
discussed elsewhere.) We see that tp = 2 for g = 3/4. 
This then establishes Jl7],[l(| that the transition is first 
order for g > 3/4. 

In the case of a first-order two-step phase-transition, 
this may not be the critical point of the phase-coexistence 
line of a many-step system. In such a case, a phase dia- 
gram of the type shown in Fig. [l]b.2 is plausible, where a 
first-order line continues from the zero-density transition 
point, ending on the coexistence curve (hidden region in 
Fig. lb. 2, see figure caption). The first-order line repre- 
sents a tranition from a uniform phase of single steps to 
a uniform phase of pairs, and, at the tranition, there will 
be diverging fluctuations in the mean spacing between 
the members of the pairs (for 3 > tp > 2). In recent 
experiments on Si(113) such a phase diagram has been 
observed Q for certain azimuthal angles of the miscut di- 
rection, though the low-density line is yet to be detected. 

Let us now consider the effect of the three-body inter- 
action [|l^,^3[|T). Even though is marginal (dimen- 
sionless at d = 1), there are contributions in v% renor- 
malization, from both vi and g (i.e. U2) via the reunion 
exponent of three steps. Note that V3 does not affect U2 
or g. Taking into account the first nontrivial contribu- 
tion (Figs. 2d, e), the recursion relation for U3 = L 1 
is (d = 1) 



v 3 



(3) 



where the numerical value of the constant A > is 
not crucial. The fixed points (Fig. lc) are U3 = 
and U3 = — 3U2/A. The stability of the — fixed- 
point depends on the sign of u^. At — it* (fixed 
point A), a small U3 is irrelevant, but the new fixed- 
point C, U3 = — 3u s /A < 0, is unstable. In contrast, at 
u 2 = u u (fixed point B), with g > 0, U3 is relevant, and 
M3 = (3|ii u |)/A > is stable. The flows for any g > are 
shown in Fig. ^c. 

Figure [j]c shows that for certain combinations of g > 
0, U2, and U3, the steps in the high-temperature phase 
behave like the 7/3 = case (region Pf ) but for sufficiently 
large u-^ < there will be triplets of steps but never pairs 
(Phase P3). The transition will be first-order but with 
singularities determined by the (/-dependent fixed point 
C of Fig. Ic. We like to associate Si(lll) vicinal surfaces 
M with this fixed point and the associated flow, because 
of similar behavior. For fixed points B and D, there could 
be pairings of steps, induced by a variation in 1*2- Such 
pairings are seen ||, e.g., in Si(001). A many-step system 
will show a tricritical behavior at these fixed points, as 
e.g. in Si(113). These two fixed points B and D of Figs. 
lc,d have the same thermodynamic exponents but they 
differ in the three-step correlation functions. Moreover, if 



U2 is a temperature-like axis, then the temperature devi- 
ation is not the proper scaling variable at D. Taking this 
as a reason for the large range j| in the observed expo- 
nent /3, we may associate Si(l 13) with the fixed point D 
with g around 3/4 and Si(001) with fixed point B. There 
is, in fact, no requirement that the long-range interaction 
should be exactly at g — 3/4 for Si(113). Unfortunately 
no systematic studies are available for the dependence of 
g on the various miscut parameters. In case g depends 
on the azimuthal angle even slightly, it is possible that 
for certain azimuthal angles g > 3/4 and for these cases, 
one would see a phase diagram of type Fig. lb. 2, while, 
for cases with g < 3/4, one would see type Fig. lb.l. 

We now make a conjecture that all vicinal surfaces are 
described by the RG flows and the fixed lines we ob- 
tained in the three-dimensional parameter space. Other 
higher-order short-range interactions are irrelevant. All 
universal properties are determined by the fixed lines, 
and it is the long-range parameter g that determines the 
universality class. 

Let us now go back to Eq. (|l|) as a quantum- 
mechanical problem. With V3 = 0, the connection of 
this Hamiltonian with the Calogero-Sutherland model 
[ ]l4j has been noted in the past |l5|,|| . The quantum prob- 
lem is described by the stable fixed line and the negative- 
g-part of the unstable line in Fig. la. Furthermore, the 
Calogero-Sutherland model describes ^] a gas of nonin- 
teracting particles obeying Haldane-Wu fractional exclu- 
sion statistics 23 1 . The long-range 1/r 2 potential is best 



thought of as a statistical interaction, with g determin- 
ing the statistics of the particles. We have established 
this directly in our approach by a computation of the 
quantum second virial coefficient along the fixed lines 
of Fig. la, since the second virial coefficient was shown, 
from general considerations, by Murthy and Shankar |2^] 
to determine the statistics uniquely. We wish to dis- 
cuss this technical issue elsewhere. Suffice it to say here 
that the non-renormalization of g in our RG approach 
is a direct manifestation of the fact that the statistics 
of the particles is independent of length scale and hence 
a renormalization- group invariant. This correspondence 
therefore tells us that the steps on vicinal surfaces are an 
analog realization of one-dimensional fractional-statistics 
particles. Though fcrmions are extensively used to study 
the equilibrium properties of steps, the latter are bet- 
ter represented by fractional-statistics particles with only 
short-range interactions. 

Our main results have already been summarized. It 
would be interesting to study the case of unequal elastic 
constants of the steps, and the regions not describable by 
any fixed point, where other details of the surface, like 
lattice periodicity, may be important [ jTl| . We end with 
the suggestions that (i) attempts be made to determine 
accurately the long-range (or statistics) parameter g and 
its dependence on the miscut features, as e.g. azimuthal 
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angle, for various surfaces and materials, (ii) measure- 
ments be done for three-step correlations, and (in) the 
thermodynamic behavior of a finite number of steps be 
studied even if in finite geometries (by putting barriers). 
Most interesting would be the special case of g = — 1/4 
which, in the surface context, is the point of coalescence 
of the stable uniform-density-phase fixed point and the 
unstable critical point, and, in the quantum context, is 
the semion, exactly half way between boson and fermion. 
Can one find such a surface of a metal? 
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bliauacliarjec-niukheiji figl 

FIG. 1. (a) Ronormalization-group fixed-lines for V3 = 0. These are stable (thick line), unstable-critical (dotted), and 
unstable-first-order (hatched), (b) A schematic orientational phase diagram in the T vs p plane for 112 < 0. For g < 3/4, 
there is a tricritical point (b.l) with /3 ~ .5. For g > 3/4, a low-density first-order line connects the p = transition to the 
coexistence curve either at the p 7^ critical-point of the latter or at an off-critical point (b.2). (c) The fixed points in the 
(w2,M3) plane for a g > 0. Arrows indicate the RG flows. The phases are Pf (uniformly-stepped, fractional-statistics like), P3 
(a phase with 3-steps bunching but no pairs), P2 (pairs of steps but no triplets), and Pt (all possible bunches), (d) The fixed 
lines on a «2 - «3 plane as g(> —1/4) is varied. Si surfaces showing similar behavior are indicated (not necessarily having 
same g). Different types of lines represent different behaviors as noted in (a) and (c). ( Fig. 3 shows the fixed lines in a 
three-dimensional (g, 112, 113) plot.) 
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FIG. 2. (a) A typical configuration of two steps consisting of sequences of bound (A) and unbound (B) fragments of various 
lengths, (b, c) A zeroth-order and a first-order reunion diagrams of two steps, (d) and (e) show the first- and second-order 
diagrams for V3 renormlization. 
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